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Abstract. We present the soft X-ray spectrum of the LINER galaxy M81 derived from a long observation with the 
XMM-Newton RGS. The spectrum is dominated by continuum emission from the active nucleus, but also contains 
emission lines from Fe L, and H-like and He-like N, O, and Ne. The emission lines are significantly broader than 
the RGS point-source spectral resolution; in the cross dispersion direction the emission lines are detected adjacent 
to, as well as coincident with, the active nucleus. This implies that they originate in a region of a few arc-minutes 
spatial extent (1 arc-minute ~ 1 kpc in M81). The flux ratios of the OVII triplet suggest that collisional processes 
are responsible for the line emission. A good fit to the whole RGS spectrum is obtained using a model consisting 
of an absorbed power law from the active nucleus and a 3 temperature optically thin thermal plasma. Two of the 
thermal plasma components have temperatures of 0.18 ± 0.04 keV and 0.64 ± 0.04 keV, characteristic of the hot 
interstellar medium produced by supernovae; the combined luminosity of the plasma at these two temperatures 
accounts for all the unresolved bulge X-ray emission seen in the Chandra observation by Tennant et al. (20011. 
The third component has a higher temperature (1.71q 5 heV), and we argue that this, along with some of the 0.64 
keV emission, comes from X-ray binaries in the bulge of M81. 

Key words. X-rays: galaxies - ISM: supernova remnants - Galaxies: individual: M81 - Galaxies: active 



1. Introduction 

M81 is an Sab spiral galaxy hosting a low luminos- 
ity Seyfert nucleus which shows the characteristics of 
a "low ionization nuclear emission-line region" (LINER, 
Hcckman ll980fl . LINERs make up a significant fraction of 
all galaxies (between 1/5 and 1/3; Ho et al. H997[) . and be- 
cause M81 is the nearest LINER it is an important target 
with which to investigate the X-ray emission from such 
objects. 

M81 has been the subject of a number of X-ray stud- 
ies. It was first observed with Einstein (Elvis & van 
Speybroeck 1982, Fabbiano 1988) which resolved several 
discrete sources in M81, the brightest of which was associ- 
ated with the active nucleus. An apparently diffuse emis- 
sion component associated with the bulge of M81, spa- 
tially extended over a few arc-minutes was detected first 
with ROSAT (Roberts & Warwick 12771701 Immler & Wang 
2001} and confirmed using Chandra data (Tennant et al. 
I2001f) . Low resolution X-ray spectroscopy of M81 from 
the Einstein MPC and IPC suggested that the nuclear X- 



ray source has an intrinsic column density of > 10 21 cm 2 
(Elvis & van Speybroeck Fabbiano H555|l . Later, 

spectra from ROSAT (Immler & Wang HHl, BBXRT 
(Petre et al. HUMjl . ASCA (Ishisaki et al. IT9UrJjl and 
BeppoSAX (Pellegrini et al. I2l70"0"jl indicated the presence 
of a soft (< 1 keV) thermal component in addition to an 
absorbed power law component (r ~ 2) from the nucleus. 
In this paper we present the soft X-ray spectrum of M81 
at much higher resolution, from an observation with the 
XMM-Newton reflection grating spectrometer (RGS, den- 
Herder et al. I 
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* Based on observations obtained with XMM-Newton, an 
ESA science mission with instruments and contributions di- 
rectly funded by ESA Member States and the USA (NASA) 



2. RGS data reduction and spectral analysis 

2.1. Data reduction 

M81 was observed by XMM-Newton on the 22nd and 23rd 
April 2001 for a total of 138 ks. The RGS data were 
reduced using the XMM-Newton science analysis system 
(SAS) public release version 5.2 and the latest calibration 
files (as of December 2001). The source spectrum was ex- 
tracted from a region centred on the nucleus, enclosing 
90% of the point spread function in the cross dispersion 
direction, while first and second order selection was per- 
formed so as to include 93% of the expected CCD pulse 
height distribution of the source photons. Background re- 




gions were selected in the cross-dispersion direction so as 
to exclude 99% of the nuclear point spread function; iden- 
tical order selection was performed on nuclear and back- 
ground regions. Instrumental features such as that near 
the oxygen edge at ~23 A were corrected by dividing the 
effective area calibration in each response matrix by the 
model/data ratio from a power law fit to the RGS spec- 
trum of the continuum source Mrk 421 (this correction 
is < 10% over most of the wavelength range). First and 
second order spectra from RGS1 and RGS2 and their cor- 
responding response matrices were resampled to match 
the RGS1 first order spectrum, combined, and rebinned 
by a factor of 3 (to ~ 30 mA bins) to improve signal to 
noise ratio before spectral analysis using XSPEC. 

2.2. X-ray emission line gas 

Fig- \H shows the RGS spectrum. A number of prominent 
emission lines from H-like and He-like N, O, and Ne as 
well as L shell lines from Fe XVII - Fe XXIII are visible 
above the continuum; these are labelled in Fig.^ and the 
most significant lines are listed in Table. The lines are 
considerably broader than the RGS line spread function 



Table 1. Fluxes and equivalent widths of the most signif- 
icant emission lines in the RGS spectrum. The fluxes have 
been corrected for Galactic absorption. Errors are quoted 
at 95% for 1 interesting parameter. 
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Ne IX" 
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Fe XVII 
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Flux 



(10~ 5 photon cm" 2 s" 1 ) 
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-0.8 
1 1+ 07 



2.0Z 



-55 
40+48 
^-34 



0.53 
0.83 
1.4 
2.5 
1.3 
2.0 
1.1 



o.s 

+0.41 
0.47 
+0.40 
0.56 
+0.6 
0.7 
+0.6 
0.7 
1.3 
1.2 
1.4 
1.4 
1.1 
0.8 



* blended with emission from Fe XXIII 

** blended with emission from Fe XIX- XXI 

for a point source (~ 70 mA FWHM; den Herder et al. 
2001). For example, O VIII Lya at ~ 19A is expected 
to have relatively little contamination from neighbouring 



M.J. Page et al.: X-ray emission line gas in the LINER galaxy M81 



3 



lines, but a gaussian fit to this line shows that its intrinsic 
width is inconsistent with a = at > 99%. This level of 
line broadening could occur if the gas is extended over a 
region of ~ 2 arc-minutes, because the RGS is a slitless 
dispersive spectrograph and its resolution is degraded for 
extended objects. Alternatively, the broadening could be 
due to kinematic motions in the gas, for example if it is 
associated with the active nucleus. 

To distinguish between these alternatives, we have ex- 
amined the spatial distribution of the emission line gas in 
the RGS cross dispersion direction by constructing spectra 
in 40 arc-second wide strips at a range of distances from 
the nucleus of M81 (for comparison, the cross dispersion 
region used in Sect. 12. l1 varies from ~ 52 arc-seconds wide 
around 20A to ~ 75 arc-seconds wide near the ends of 
the spectrum). The spectra were binned to 60 mA and a 
background spectrum was constructed from the average 
of the two most off-axis strips; this was subtracted from 
all the spectra. The emission-line rich 14-20 A parts of 
these spectra are shown in Fig. [21 The emission line gas, 
and in particular O VIII Lya, is extended over more than 
an arc-minute implying that the excess broadening of the 
emission lines is due to their spatial extent, rather than 
kinematic motions. 

We further investigate the nature of the emission 
line gas using plasma diagnositics from the He-like OVII 
triplet, as described in Porquet & Dubau (2000). Taking 
the 21-23 A portion of the RGS spectrum, and assum- 
ing an underlying power law continuum (r = 2, see Sect. 
I2.3fl . we fitted the OVII lines as 3 gaussians with the rel- 
ative wavelengths of the resonance (w), intercombination 
(x+y) and forbidden (z) lines fixed at their theoretical ra- 
tios. The best fit has an acceptable x 2 / u — 56/51, and 
is shown superimposed on the 21-23 A RGS spectrum in 
Fig. UJl. In Fig. [SJd we show the confidence contour of the 
(x+y+z) against w line strength, a diagnostic of the ion- 
ization mechanism. The solid line shows the line ratio G = 
(x+y+z)/w = 4. As shown in figure 7 of Porquet & Dubau 
(2000), plasmas in which photoionization is the dominant 
ionization mechanism are expected to have ratios to the 
left of this line (i.e. G > 4); this is ruled out at 95% confi- 
dence (A% 2 > 4 for one interesting parameter G) in M81. 
For plasmas with line ratios to the right of this line, such 
as that observed in M81, collisional excitation is impor- 
tant and may be the dominant emission mechanism. If the 
O VII lines arise in a collisionally ionized plasma, the G 
ratio provides a temperature diagnostic; by comparison 
with Porquet et al. <|2001[1 we find that the lines arise in a 
plasma with kT < 0.26 keV at 90% confidence. 

Fig. |3t shows confidence contours of the density diag- 
nostic z against (x+y). The line shown on the plot corre- 
sponding to ratio R = z/(x+y) = 1 is excluded at > 95% 
along with all of parameter space to the right of this line, 
implying (see Porquet & Dubau 2000 figure 8) that the 
line emitting plasma has a density of n e < 10 1:L cm -3 . The 
best fit has very weak x+y lines relative to z, and conse- 
quently lies to the left of the line R = 3.5, hence the line 
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Fig. 2. Spectral cuts labeled with distance from the nu- 
cleus in the cross dispersion direction (XDSP) demon- 
strating that the line emission comes from a spatially ex- 
tended region. The dashed lines correspond to prominent 
emission lines from FeXVII and OVIII. 



ratio observed in M81 is best reproduced by a low density 
(n e < 10 9 cm -3 ) plasma. 

2.3. Modelling the whole RGS spectrum 

It is evident from Fig. ^ that the emission lines in the 
soft X-ray spectrum of M81 are superimposed on a strong 
continuum source, almost certainly from the active nu- 
cleus at the heart of M81. Therefore we began by fitting 
a power law model with absorption (using the XSPEC 
model TBabs of Wilms et al. ESMij) by the line of sight 
Galactic column density Nh = 4.16 xlO 20 cm~ 2 (Dickey 
& Lockman[1990). The parameters of this fit, and subse- 
quent fits, are given in Table |3 We found that the power 
law model is a poor fit with x 2 l v = 1289/993; the model 
overpredicts the data at both the high and low wavelength 
ends of the spectrum, and has too small an O I edge at 
23 A. We therefore added additional cold absorption, and 
obtained a much better fit (x 2 jv = 1085/992), with an 
additional Nh ~ 5 xlO 20 cm~ 2 . This additional column 
density confirms, and is consistent with, previous results 
(e.g. Pellegrini et al. 2000) that the nucleus of M81 has 
some intrinsic absorption. 
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Fig. 3. Panel a. shows the O VII triplet and best fit 
model. Panel b. shows the 68% (solid), 90% (dashed) 
and 95% (dotted) confidence contours for the intercom- 
bination (x+y) and forbidden (z) line strength against 
the resonance (w) line strength; photoionization powered 
lines should have line ratios to the left of the line G = 
(x+y+z)/w = 4. Panel c. shows 68% (solid), 90% (dashed) 
and 95% (dotted) confidence contours for the forbidden 
(z) line strength against the intercombination (x+y) line 
strength; the ratios R = z/(x+y) = 1 and R = 3.5 corre- 



spond to densities of < 
tively. 



10 9 cm 3 respec- 



We then attempted to model the emission line com- 
ponent of the spectrum. As shown in Sect. 12.21 the O VII 
triplet is consistent with an origin in a collisionally ionized 
gas, and the line emission comes from a region extended 
over several kpc. We therefore began modelling the emis- 
sion lines by adding a single-temperature, solar-abundance 
'Mekal' thermal plasma model to the absorbed power law 
model. This resulted in a significant improvement in the fit 
( x 2 /v = 876/990) for a plasma temperature of 0.66 + 0.04 
keV. 

However, the emission lines are too narrow in this 
model (the excessive width of the emission lines has al- 
ready been pointed out in Sect. l2~2l . and some emission 
lines, particularly those from OVII, have much lower in- 
tensities in the model than in the observed spectrum (see 
Fig. We addressed the first of these two problems by 
adding 2 extra Mekal components, shifted by ±0.5% in 
wavelength, with temperatures tied to that of the first 
Mekal component, but with normalisations allowed to 



vary. This broadened the line profile of the model to mimic 
the degraded resolution of the RGS for the extended emis- 
sion line region. We then added a second, lower temper- 
ature, broadened Mekal component to the model. This 
resulted in another significant improvement in x 2 , repro- 
ducing most of the emission lines well, with the second 
plasma component at a best fit temperature of ~ 0.2 keV, 
consistent with the temperature deduced from the O VII 
lines (Sect. However, although x 2 is good, the blend 
of Ne X with Fe XXI-XXIII at - 12.1 A is still underpro- 
duced by this model, and so a third, higher temperature, 
broadened Mekal plasma was added. This resulted in yet 
another significant improvement in the fit, to a final best 
X 2 /v = 803/980, with a best fit temperature for the third 
Mekal component of 1.7^5 5 keV. The model reproduces 
the spectrum well, as shown in FigQ] 

We have not attempted to derive elemental abun- 
dances for the soft X-ray emission line gas, because the 
abundance ratios are somewhat degenerate with the tem- 
perature structure of the gas. However, we have verified 
that our assumption of solar abundance is reasonable as 
follows. Fe and O are responsible for the strongest soft 
X-ray lines in M81, so the ratio of Fe to O will be the 
most important abundance ratio in terms of the resultant 
emission line spectrum, and the easiest to determine. We 
have therefore repeated the three-Mekal fit, once with the 
Fe abundance halved, and once with the Fe abundance 
doubled, but with the abundances of other elements (in- 
cluding O) fixed at the solar values (Anders & Grevesse 
1989}. We find that with the Fe abundance halved, the 
best fit is worse by A% 2 = 4, while doubling the Fe abun- 
dance results in a best fit that is worse by A% 2 = 9, with 
respect to the three-Mekal, solar-abundance fit in Tabled 
This suggests that the thermal plasma in M81 has an Fe 
to O ratio reasonably close to solar. 

3. Discussion 

In Sect. 12.21 we showed that the soft X-ray emission lines 
in M81 come from a region which is extended in both 
the dispersion direction (evidenced by the excess width of 
the lines) and cross dispersion direction (see FigQJ of the 
RGS. We also showed that the diagnostic OVII triplet is 
emitted by a low density plasma and is not powered by 
photoionization from the central AGN. Furthermore, in 
the previous section we showed that a 3 temperature ther- 
mal plasma model provides a good description of the emis- 
sion line component of the M81 RGS spectrum in conjunc- 
tion with absorbed power law emission from the nucleus. 
Two of the plasma components have well constrained tem- 
peratures of 0.18 + 0.04 keV and 0.64+0.04 keV, which are 
characteristic of the hot interstellar medium produced by 
supernova explosions; in this case we would expect these 
emission lines to come from a genuinely diffuse region (or 
regions) in the bulge of M81. 

X-ray imaging studies, first using ROSAT (Roberts 
& Warwick 123001 Immler & WangEUffl) and later using 
Chandra (Tennant et al. 1200 1[1 have indeed found appar- 
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Fig. 4. Power law x intrinsic neutral absorption + 3 Mekal model fit to the RGS spectrum. The top panel shows the 
data (grey points) and model (bold black line) while the bottom panel shows the data/model ratio. Also shown in 
the top panel is the contribution from the three different temperature thermal plasma components: 0.18 keV (thin 
solid line), 0.64 keV (dashed line) and 1.7 keV (dotted line). Inset a displays data and single-temperature Mekal 
model around the OVII triplet, showing that O VII lines are not reproduced by the single temperature model. Inset 
b shows data and two-temperature Mekal model around the Ne X / Fe XXI- XXIII blend, indicating that a third, 
higher-temperature model component is required to reproduce these lines. 



ently diffuse emission associated with the bulge of M81. 
Furthermore, Immler & Wang (2001) extracted a PSPC 
spectrum of the M81 bulge using a 1 arc-minute radius 
source circle, including both the active nuclues and the 
diffuse emission. They found a good fit using an absorbed 
power law to represent the nucleus and a two tempera- 
ture thermal plasma to represent the apparently diffuse 
emission. Their best fit plasma temperatures, 0.15 ± 0.02 
keV and 0.63±0.11 keV, are in extremely good agreement 
with the temperatures obtained from the RGS spectrum. 
However the spatial resolution of ROSAT (~ 5 arc-second 
for the HRI and ~ 20 arc-second for the PSPC) was such 
that the emission from a large population of weak point 
sources could easily be confused with genuinely diffuse 



emission. 



The much higher resolution observation performed 
with Chandra provides a much more accurate determi- 
nation of the truly diffuse emission from the nucleus. 
Tennant et al. I|2001fl found that of the bulge emission, 
even excluding the nucleus, 64% comes from the resolved 
point source population. The remaining unresolved bulge 
emission reported by Tennant et al. (|2001|l has a coun- 
trate of 0.092 s" 1 in the Chandra ACIS-S3 chip. The 0.18 
keV and 0.64 keV thermal plasma components from our 



best fit model in Table 1 would correspond to countrates 
of 0.036 s _1 and 0.158 s _1 in the same chip. The lowest 
temperature component alone is unable to account for all 
the unresolved bulge emission, but the combinaton of the 
two exceeds it. 

This implies that some of the emission associated with 
the 0.64 keV thermal plasma, and all of the emisson asso- 
ciated with the highest of the three plasma temperatures 
(l-7ip'5 keV) must arise in part of the bulge population re- 
solved by Chandra. We therefore propose low-mass X-ray 
binaries as the most likely origin of this higher tempera- 
ture emission. 

The line-of-sight velocity dispersion in the bulge of 
M81 is 165 ± 10 km s" 1 (Prichet IT97gl Pellet & Simien 
1982$, so the 0.2 keV plasma has similar energy per unit 
mass to the stellar component and is bound to the bulge. 
Using the cooling curves from Landini & Monsignori Fossi 
( 1990) we find that the cooling time for the 0.2 keV com- 
ponent is ~ 4 x 10 7 / y/f years where / is the filling factor 
of the gas. This is much shorter than the time since the 
last M81-M82 perigalactic passage (~ 5 x 10 s years ago, 
Brouillet et al. 119911 Chandar et al. 12001(1 and hence the 
0.2 keV gas, which has a mass of W 6 /VJ Mq, cannot be a 
relic of a nuclear starburst induced by the passage of M82 
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Table 2. Spectral fitting of the whole M81 RGS spectrum. 
All models include Galactic absorption of Nh = 4.16x 10 20 
cm -2 



(1) 

Model 


(2) 
T or kT 


(3) 
flux 


(4) 
N H 


(5) 

2 / 

X /v 


PL 


1.67 ±0.02 


105 ±1 




1289/993 


PLxNh 


2.06 ±0.06 


99.9 ±3.8 


4.8 ±0.8 


1085/992 


PLxNh 
+Mek 


1.94 ±0.07 
0.66 ± 0.04 


95.8 ±3.7 
4.31°;? 


3.3 ±0.8 


876/990 


PLxNh 
+Mek3 


1.92 ±0.07 
0.64 ±0.04 


95.4 ±3.7 

4 7 +o.9 


3.1 ±0.8 


866/988 


PLxNh 
+Mek3 
+Mek3 


1.93 ±0.07 
20+ 04 
0.66 ± 0.05 


94.5 ±3.8 

1.5 ±0.5 
4 9+0. s 


3.4 ±0.8 


821/984 


PLxNh 
+Mek3 
+Mek3 
+Mek3 


1.94 ±0.06 
0.18 ±0.04 
0.64 ±0.04 
1 7 +2A 

*" 1 -0.5 


91.6±H 

1 1+0.5 
1 - l -0.4 

3-4t^ 


3.4 ±0.8 


803/980 



Explanation of columns: 

1 PL = power law, 

Nh = neutral absorber ('tbabs' model in XSPEC), 
Mek = thermal plasma ('Mekal' model in XSPEC), 
Mek3 = thermal plasma ('Mekal' model in XSPEC) with 
components shifted by ±0.5% in wavelength to mimic the 
broadening of the lines produced by the spatial extent of 
the emission region. 

2 power law photon index or temperature of thermal 
plasma (keV) 

3 flux of model component in the energy range 0.3-2.0 keV, 
in units of 10 -13 erg cm -2 s -1 without Galactic 
absorption of N H = 4.16 x 10 20 cm" 2 

4 intrinsic neutral column density (10 20 cm -2 ). 



because the gas must have been replenished or reheated 
much more recently. The 0.64 keV component has a mass 
of 1.4 x lO 6 /-^// Mq, and potentially has a much longer 
cooling time, ~3x lO 5 /-^/ years, but it is too hot to be 
bound by the virial mass of the bulge. The sound crossing 
time for the bulge of M81 is less than 10 7 years, so the 
0.64 keV gas must also be replenished or reheated on a 
timescale of a few xlO 7 years. 

The simulations of Shelton ( 193SJ) allow us to estimate 
the approximate supernova rate that would sustain the 
diffuse X-ray luminosity of the bulge of M81. The com- 
bined flux from our two lowest temperature thermal com- 
ponents, compared to the predicted ROSAT countrates 
from section 7.4 of Shelton l|1998[l corresponds to super- 
nova rates of between ~4x 10 -3 year -1 and ~4x 10 -2 
year -1 , where the lower estimate comes from the predicted 
countrate in the lower energy Rl ROSAT band and the 
higher estimate comes from the predicted countrate in the 



higher energy R4 ROSAT band. The flux in the ROSAT 
R4 band is dominated by the 0.64 keV plasma; if we as- 
sume that half of this emission is due to resolved sources 
then the R4-band estimate of the supernova rate is re- 
duced to ~ 2 x 10 -2 year -1 . If the supernova rate in the 
bulge of M81 is ~ 4 x 10 -3 year -1 or higher we expect 
there to be tens of supernova remnants within the bulge 
with ages < 10 4 years, although in the gas-poor environ- 
ment of the bulge they may have very low X-ray lumi- 
nosities. There are 41 sources detected by Chandra within 
the bulge of M81, of which ~ 30% have spectra which are 
softer than typical X-ray binaries (Tcnnant ct al. 2001). 
The collective flux of these sources is sufficient to con- 
tribute significantly to the 0.64 keV plasma emission; fur- 
ther Chandra observations will be required to determine 
whether these sources are persistent and thus potential 
supernova remnants. 

UV imaging carried out with the Hubble Space 
Telescope failed to reveal any individual massive stars in 
the bulge (Devereux et al. H997J> . and this lack of a young, 
massive stellar population leaves only type la supernovae 
as the possible origin of the X-ray emitting gas. We ob- 
tain the canonical 'expected' rate for type la supernovae 
in M81 from table 8 of Van den Bergh & Tammann <|1991[l . 
Assuming a distance of 3.6 Mpc, B magnitude of 7.31 after 
correction for Galactic extinction (de Vaucouleurs et al. 
ITMl"|) and H = 75 km s -1 Mpc -1 , we expect - 6 x 10 -3 
type la supernovae year -1 for the whole galaxy, consis- 
tent with the rate of supernovae required to support the 
diffuse X-ray emission. 

The diffuse X-ray emission may well relate to the opti- 
cal line emission which gives M81 its LINER classification. 
The bulge of M81 is bright in Ha emission (Devereux et al. 
1995] Greenwalt et al. H998[) . and in the absence of young 
massive stars Devereux et al. ©97) conclude that the dif- 
fuse Ha emission must be powered either by old post-AGB 
stars or by shocks, both of which are compatible with the 
observed UV surface brightness. If the diffuse X-ray emis- 
sion is indeed produced by supernova remnants, then fast 
shocks must be propagating through the bulge. 

There are three pieces of evidence that suggest shocks 
related to the X-ray emitting gas, rather than post-AGB 
stars, ionize the diffuse optical emission line gas. Firstly, 
the Ha emission region has a similar spatial extent to the 
X-ray line emitting plasma, but it has an asymmetric, pos- 
sibly spiral structure (Devereux et al. [T995| Greenwalt et 
al. I1998|l . and hence does not trace the underlying distri- 
bution of stars in the bulge, as might be expected if the 
gas is ionized by old post-AGB stars. Secondly, the ion- 
ized Ha emitting gas shows complex, non-rotational mo- 
tion of up to 200 km s -1 , as well as rotation at up to 300 
km s -1 (Goad lrl)76fl . The non-rotational motion could be 
driven by the expansion of hot bubbles and suggests me- 
chanical heating of the gas. Thirdly, the optical emission 
line ratios suggest that photoionization by hot stars is 
not the ionization mechanism in the bulge. For example, 
the [SII]/(Ha±[NII]) ratio is much higher in the bulge 
of M81 than in any of the HII regions, (Greenawalt et 
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al. I1998[) . suggesting that photoionization by hot stars is 
not the ionization mechanism in the bulge. However shock 
heating, with shock velocities of a few hundred km s _1 , 
can result in large [SII]/Ha ratios (Baum et al. ] 1992 
Dopita & Sutherland ll995|) as observed. Furthermore, the 
strength of [O III] A5007 observed in the central few hun- 
dred pc of M81 cannot be reproduced by photoinization 
from hot stars (Golev et al. 119961 Wang et al. I1997|) . but 
the [OIII]/H/3 ratio is similar to the predictions of the 
shock model of Dopita & Sutherland (1995 ) for shock ve- 
locities of 300 km s -1 . In this model, [OIII] A5007 is emit- 
ted by gas which is photoionized by the shock-heated ma- 
terial. For a shock velocity of 300 km s _1 the shocked ma- 
terial has an electron temperature similar to the 0.2 keV 
plasma found in the RGS spectrum (Dopita & Sutherland 
1996). This suggests that the optical line-emission is re- 
lated to the soft X-ray line-emitting plasma, and that the 
ISM in the bulge of M81 is shock-heated, possibly by su- 
pernovae. 

4. Conclusions 

The X MM- Newton RGS soft X-ray spectrum of M81 
shows emission lines from H-like and He-like N, O and 
Ne as well as Fe-L lines superimposed on a strong contin- 
uum from the nucleus. The excessive width of the emission 
lines in the dispersion direction, and their detection out- 
side the nucleus in the cross dispersion direction implies 
that the emission lines originate in a region of a few arc- 
minutes spatial extent, corresponding to a few kpc in M81. 
The OVII triplet line ratios suggest that photoionization 
is not the main ionization mechanism and that collisional 
processes must be important in producing the observed 
lines. The RGS spectrum can be fitted with a model con- 
sisting of an absorbed power law from the nucleus and 
a 3 temperature optically thin thermal plasma. Two of 
the thermal components have temperatures (0.18 ± 0.04 
keV and 0.64 ± 0.04 keV) which are consistent with the 
hot interstellar medium produced by supernovae; the com- 
bined flux from these two components fully accounts for 
(in fact exceeds) the unresolved bulge emission seen with 
Chandra. The X-ray emission could be produced by type 
la supernova rates of ~ 4 — 20 x 10 -3 year -1 , which is not 
unreasonable for M81. We propose that the shocks gener- 
ated by the supernova remnants could also be responsible 
for the observed optical line emission in the bulge of M81. 
The third X-ray emitting thermal plasma component has 
a higher temperature (1-7_q\ H keV) and in order not to ex- 
ceed the unresolved X-ray emission found with Chandra 
we propose X-ray binaries in the bulge of M81 as a likely 
origin of this emission, as well as ~half of the 0.64 keV 
emission. 
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